EeV neutrinos associated with UHECR sources 
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Electromagnetic energy losses of charged pions and muons suppress the expected high energy, 
> lO'^^ eV, neutrino emission from sources of ultrahigh energy, > 10^^ eV, cosmic-rays. We show 
here that > lO'^^ eV photons produced in such sources by neutral pion decay may escape the sources, 
thanks to the Klein-Nishina suppression of the pair production cross section, and produce n pairs 
in interactions with the cosmic microwave background. The flux of muon decay neutrinos, which are 
expected to be associated in time and direction with the electromagnetic emission from the sources, 
may reach a few percent of the Waxman-Bahcall bound. Their detection may allow one to directly 
identify the sources of > 10^^ eV cosmic-rays, and will provide the most stringent constraints on 
quantum-gravity-induced Lorentz violation. 

PACS numbers: 95.85.Ry, 14.60.Pq, 98.70Rz, 98.70.Sa 



I. INTRODUCTION 

The observed characteristics of cosmic rays with ener- 
gies exceeding 10^^ eV [i|,0,[3| suggest that there are cos- 
mologically distributed astrophysical sources that accel- 
erate protons to this high energy [J] . The main candidate 
sources are gamma-ray bursts (GRBs) and active galactic 
nuclei ( AGN) [4] . The sources of the ultra- high energy 
cosmic rays (UHECRs) are expected to emit also high en- 
ergy neutrinos, through the decay of charged pions pro- 
duced in interactions of high energy protons with the in- 
tense radiation fields within the prospective sources @. 
Large neutrino telescopes, both operating and under con- 
struction, are expected to reach in the coming few years 
the sensitivity required for the detection of the predicted 
high energy neutrino signal [6| . 

The detection of high energy neutrinos associated spa- 
tially and/or temporally with electromagnetic emission 
from some sources will provide evidence for cosmic ray 
production in these sources. However, electromagnetic 
and adiabatic energy losses of charged pions and muons 
are expected to suppress the neutrino emission at high 
energy, where the charged particle life time exceeds its 
energy loss time [3,|Ml- For example, in the case of GRBs 
the neutrino flux is expected to be suppressed above a 
few PeV [1, Q to a level which would make it difficult 
to detect. This would make it difficult to obtain direct 
evidence for acceleration of protons to energy > 10"'^^ eV. 

We discuss here a new type of high energy, > 10^® eV, 
neutrino signal, which is expected to be associated with 
UHECR sources, and which is not suppressed by the en- 
ergy losses of charged pions and muons. Neutral pions 
produced in p7 interactions decay to photons before un- 
dergoing significant adiabatic and/or electromagnetic en- 
ergy loss, since their life time is much shorter than that 
of charged pions and since they are not trapped by mag- 
netic fields and hardly interact with the plasma and with 
the radiation field. The resulting high energy, ~ 10^^ eV, 
photons are likely to escape the source, due to the Klein- 
Nishina s upp ression of the pair production cross section 
(see e.g. [1^ for a discussion of the GRB case), and may 



then produce muon pairs in interactions with cosmic mi- 
crowave background (CMB) photons. Muon decay leads 
to the production of high energy, > 10^^ eV, neutrinos, 
which are associated, both in direction and in time, with 
the electromagnetic emission from the sources. High en- 
ergy neutrino telescopes like ARIANNA 11] may allow 
one to detect these > 10^^ eV neutrinos. 

In what follows we first discuss fi pair production in 
interactions of > 10^^ cV photons with CMB photons, 
and show that ~ 5% of the photon energy is expected to 
be converted to (prompt) high energy muons. We then 
discuss the production and escape of > 10^^ eV pho- 
tons from GRB sources, as an example of likely sources 
of UHECRs which despite being optically thick to pair 
production for 1 GeV to 1 PeV photons, are likely to 
be optically thin for > 10^^ eV photons. Finally, we dis- 
cuss the expected neutrino flux and spectrum for UHECR 
sources which, like GRBs, deposit via pion production a 
significant fraction of the energy of UHECR protons in 
high energy photons, and which are optically thin for 
> 10^^ eV photons. 



II. fi^ PAIR PRODUCTION IN INTERACTIONS 
WITH CMB PHOTONS 



Consider a photon of energy e propagating through an 
isotropic radiation field with photon number density per 
unit energy dn/de. The (inverse) mean free path (mfp) 
for e^ or /i pair production is [8| 
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Here a{E) is the cross section for interaction with en- 
ergy E in the center of momentum frame, and m = rrie 
{m^) for e (/i^) production. For a thermal radiation of 
temperature T, dn/de = 7r~^(ftc)~'^e^/[exp(e/r) — 1], we 



have, using u — E^/ATe, 



ri = 



2^3 



TT^ihc) 



'1-e- 



du ua[\/Au£T) In 

'(mc2)2/Te 

(2) 
At low energy, e < (rnc^)^/3rcMB ~ 10^^ eV, muon 
pair production is possible only in interaction with pho- 
tons from the high energy tail of the CMB photon dis- 
tribution, and the e* pair production mfp, Ze, is much 
shorter than the jjl pair production mfp, l^. The two 
mfp's become comparable at higher energy (note that the 
cross sections for /i* and e production are similar for 
E > 2m,,c2). Using eq. (0) we find R = 1^,/le = 36, 12, 8.6 
for Te/{m^c^Y — 1:3,5 respectively, corresponding to 
£ = (0.4, 1, 2) X 1020(1 -f z)-i eV (Using T = rcMB(^) = 
2.3 X 10^'*(1 + z) eV). For photons in the energy range 
of e « 10i''(l + z)~i eV to e « 10^0(1 + z)-^ eV, the 
fraction of high energy photons converted to /i pairs 
is thus ff, = 1/(1 -t- i?) « 0.03 to 0.1, and the mfp is 
l^Tn 10^^(1 -I- z)-3 cm. 

The fraction of photon energy converted to muons is 
higher than /^, due to the following reason. The cen- 
ter of momentum energy of the pair production interac- 
tion is well above the electron rest mass. This implies 
that the e produced are relativistic in the center of 
momentum frame, and that most of the initial photon 
energy is carried (in the CMB frame) by one of parti- 
cles produced. The resulting high energy electron (or 
positron) would lose most of its energy in a single inverse- 
Compton scattering of a CMB photon (which is deep in 
the Klein-Nishina regime), leading to the production of 
a secondary high energy photon, with energy which is 
not much smaller than that of the original photon. This 
process increases the effective mfp for photon energy loss 
to e^ to ~ lOZe (e.g. fig. 4 of [l2|)- Since the electron 
and photon mfp's are similar, this implies that muon pair 
production by secondary photons would increase the frac- 
tion of energy loss to muons by a factor of a few beyond 
/^ (The further small additional contribution to muon 
production via 67 -^ e^^ may be neglected [131) ■ ^^ 
we show below, while the prompt neutrinos (from the 
decay of prompt muons produced in an interaction of 
the original high energy photons with a CMB photon) 
are expected to arrive nearly simultaneously with the 
photons emitted by the source, neutrinos produced by 
interactions of secondary photons may suffer very long 
time delays due to deflection of the secondary high en- 
ergy electrons by inter-Galactic magnetic fields. Since 
we are interested in neutrinos which may be associated 
directly with their sources, we restrict our discussion to 
neutrinos produced by the decay of prompt muons 

For the thermal CMB spectrum, e 
in interactions with high energy, e 
is dominated by photons of energy ~ Tcmb ^ IQ^^ eV, 
well above the threshold energy, ~ (jrieC?)'^ je ~ 10"^ eV. 
However, the number density of ^ 10^^ eV background 
photons may be dominated by radio emission from ob- 
jects like galaxies and AGN, rather than by the CMB. 
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The presence of such additional "radio background" pho- 
tons may increase the e* pair production rate, and thus 
reduce the value of f^. Assuming that the ~ 10^^ cV 
(~ 10 MHz) radio background observed at Earth is of 
extra-Galactic origin, l^ is reduced, with respect to its 
value in the presence of CMB photons only, by a factor 
of ~ 2 for £ ~ IQ-'^^ eV, and by a factor of ~ 10 for 
e - 10^0 eV (see, e.g., fig. 1 of (H). The actual re- 
duction in If. is expected to be much smaller, since the 
~ 10 MHz radio background is dominated by Galactic 
emission. The extra-Galactic contribution to the radio 
background is probably lower than 10% of the observed 
background (see discussion in sec. 4 of [l^), which im- 
plies that Ze is not much affected for e < 10^° eV. More- 
over, the contribution of radio background photons to 
e pair production interactions becomes smaller with in- 
creasing redshift: At higher redshift the energy density 
of the CMB is higher, while the energy density of the ra- 
dio background is expected to be lower, since this back- 
ground is produced by sources like galaxies and AGN and 
is therefore accumulating on a time scale comparable to 
the age of the universe. Thus, for cosmic-ray sources ly- 
ing at z = 1 to z = 2, which are expected to dominate 
neutrino production (e.g. 0), we expect the effect of 
radio background photons to be small. 



III. HIGH ENERGY PHOTON EMISSION: THE 
GRB CASE 

Both GRBs and AGN may be capable of accelerating 
protons to ultra-high energy. The fact that out to ^ 
100 Mpc, the propagation distance of 10^° eV protons, 
there is no AGN luminous enough to meet the energy 
output constraint, L > lO^^erg/s for proton acceleration 
to 10^0 eV Q, suggests that AGN are not the sources of 
UHECRs observed at Earth. We therefore focus below 
on high energy photon production in GRBs. 

Observed GRB properties suggest that they may accel- 
erate protons to > 10^° eV and that they may produce 
the UHECR flux observed at Earth [H, [l^- In [1] it 
was shown that > 10^^ eV protons, which are acceler- 
ated within the gamma-ray emitting region of the GRB 
source, lose a significant fraction, ~ 20%, of their energy 
to pion production before escaping the source (on aver- 
age, a proton undergoes a single pion production inter- 
action before escaping). In fact, for protons accelerated 
to > 10^^ eV within the gamma-ray emitting region, the 
fraction of energy lost to pion production may be some- 
what higher than 20%. In 8], only the A-resonance con- 
tribution to the p7 interaction rate was considered. This 
is an excellent approximation for ^ 10"'^^ eV protons, for 
which the interaction with the dominant e ~ 1 MeV GRB 
photons is at the A-resonance. Since the spectrum of 
higher energy photons, for which the interaction is at an 
energy higher than the resonance, is steep, dn/dt oc e^~, 
the contribution of interactions with these photons to 
the p7 interaction rate is small (see discussion follow- 



ing eq. 1 of [^). For protons of much higher energy, 
e ^ 10^^ eV, the energy ca of photons for which the 
interaction is at the A-resonance is lower than 1 MeV, 
EA ^ l(10iSeV/e) MeV. Since the GRB spectrum is 
harder below 1 MeV, dn/de ex e~^, the relative contribu- 
tion of non-resonant interactions (compared to the reso- 
nant ones) is ~ (l/5)xln(e/10^^eV) (the 1/5 factor is due 
to the non-resonant cross section being ~ 5 times smaller; 
see eqs. 1 & 2 of Q)- This holds for protons of energy 
e < 10^9 eV, for which eA > 1 keV. The GRB spec- 
trum is expected to become self absorbed, dn/de oc e, at 
e < 1 keV (see below), andp7 interactions of e > 10^^ eV 
protons are dominated by the non-resonant interactions 
with > 1 keV photons. The P7 interaction rate for 
these protons is thus similar to the interaction rate of 
- IQi*' eV protons, (1/5) x ln(10i^/10iS) - 1. How- 
ever, since the fraction of energy loss per interaction is 
higher for interactions well above the resonance, the frac- 
tion of energy loss to pion production may be ~ 50% for 
£ > 10^^ eV protons. 

Let us consider next the question of whether or not 
>10^^ eV photons can escape the GRB source (see also 
|l(j|V The escape of high energy photons is limited by 
pair production interactions. The optical depth to pair 
production within the GRB source may be estimated as 
follows. It is widely accepted that GRB gamma-rays are 
produced by the dissipation of the kinetic energy of a 
highly relativistic outflow, usually referred to as a "fire- 
ball" , driven by accretion onto a compact object (presum- 
ably ^ 1 solar mass black hole) j^]. Let us denote the 
radius at which gamma-ray emission arises by i?, and the 
outflow Lorentz factor by F. The pair production optical 
depth for a photon of high energy e is given by the prod- 
uct of the pair production rate, 1/i' (e), and the dynami- 
cal time, the time required for significant expansion of the 
plasma, t'^ ~ R/Tc (primes denote quantities measured 
in the outflow rest frame), Tjj{e) ~ R/Tct'^^[e). t'^^{e) 
depends on the energy density and on the spectrum of the 
radiation. The (outflow rest frame) radiation energy den- 
sity is approximately given by U' = L/AirR^cT^, where 
L is the GRB luminosity. The GRB spectrum can typi- 
cally be described as a broken power law, dn/de oc e~^, 
with /? « — 1 at low energy, e < et, ~ 1 MeV, and 
/?« — 2ate > efc ~ 1 MeV. High energy photons 
with energy e' exceeding e^, = 2(TOeC^)^/eJ,, may produce 
pairs in interactions with photons of energy exceeding 
e' — 2{niec'^)'^/e' < ej, (the rest frame photon energy e' is 
related to the observed energy by e = Fe'). For e' < ej, we 



have dn/de' oc 



which implies that the number den- 



sity of photons with energy exceeding e' depends only 
weakly on energy. Thus, t^^ is nearly independent of 



energy for e' > e'^, t. 
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Here L — 10^^L52erg s~^, F — 10^'^F2.5, and the emis- 
sion radius R is related to the observed variability time 
Ai = 10~2Ai_2 s by At « R/2T'^c. Photons of lower 



energy, e < Sb, interact to produce pairs only with pho- 
tons of energy e' > 2(meC^)^/e' > ej,. Since the num- 
ber density of these photons drops like 1/e', we have 
T^-y(£ < Eb) ~ {e/eb)T^^{e > Sb), i.e. 
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The extension of GRB spectra to ^ 100 MeV and the 
characteristic variability time, At_2 ^ 1, imply F2.5 > 1. 
Since thermal pressure acceleration can not lead to much 
larg er Lorentz factors, F2 5 « 1 is typically adpoted (e.g. 

Eq. ([3]) implies that the pair production optical depth 
is large, ^ 10^, for high energy, £ > £5 « 10^^F| 5 eV, 
photons. However, at very high energy the source might 
become optically thin again. This is due to the fact that 
the spectra do not extend as dn/de' oc e'~^ to arbitrar- 
ily low energy. Synchrotron self-absorption is expected 
to "cutoff" the spectrum at some energy e'^ <C ej,, below 
which the spectrum is expected to follow dn/de' oc e' . 
For high energy photons with energy e' exceeding £^ = 
2(77ieC^)^/e^, pair production would be dominated by in- 
teractions with photons of energy ~ e'^, and the Klein- 
Nishina suppression of the pair production cross section 
would imply t-^^ oc 1/e'- For typical GRB parameters we 
expect El ~ 1 keV [13, [20] , consistent with BeppoSAX 
GRB spectra [21[. Using eq. ^ we thus have 
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We note here that the recent detection of (apparently) 
simultaneous optical/UV and gamma-ray emission ap- 
pears to be inconsistent with the standard fireball model, 
since the optical/UV flux is much higher than expected in 
a model where synchrotron self-absorption is important 
below ~ 1 keV. We have shown, however, that the op- 
tical/UV emission may in fact be naturally produced in 
this model by " residual" dissipation of kinetic energy at 
radii larger than that of the gamma-ray emission region 



IV. THE EXPECTED NEUTRINO SIGNAL 

Let us assume that the sources of UHECRs deposit a 
significant fraction of the UHEGR energy output into 
pion production, and that they are optically thin for 
> 10^^ eV photons. The fraction of UHECR energy that 
would be converted to high energy neutrinos, via the es- 
cape of high energy photons and their interaction with 
GMB photons, is w (1/2) x 0.05 x 2/3 - 2%. The first fac- 
tor accounts for the fraction of energy deposited in neu- 
tral (as opposed to charged pions), the second accounts 
for the fraction of photon energy deposited in (prompt) 
muon pairs in interaction with GMB photons, and the 
last term accounts for the fraction of energy deposited 
into neutrinos by the muon decay. This implies that the 



expected neutrino flux mav reach a few percent of the 
Waxman-Bahcall bound [9|, which is the neutrino flux 
that would be obtained if all the energy produced by the 
sources of UHECRs were deposited into pions. 

At redshift z, the energy of photons capable of pro- 
ducing muon pairs in interactions with CMB photons of 
energy ~ STcmb is ^ 10^^(1 + z)^^ eV. The character- 
istic energy of neutrinos arising from such interactions is 
~ 3 X 10^^(1 + z)~^ eV, and the energy with which they 
wih be observed at Earth is ~ 3 x 10"[(1 -h 2)/3]"^ eV. 
The neutrino flux from muon production on the CMB 
is expected to be strongly suppressed below this energy. 
The muon decay neutrino flavor ratio, $° : $° : <1>^ — 
1 : 1 :J), is expected to be modified by neutrino oscilla- 
tions [221 to $^^ : $1.^, : ^^^ = 0.4 : 0.3 : 0.3 at Earth, 
slightly different from the 1:1:1 flavor ratio for neu- 
trinos produced by charged pion decay. The predicted 
neutrino flux may be detectable by the suggested AR- 
lANNA neutrino telescopes [ll|, which is designed for 
detection of neutrinos of energy > 10^^'^ eV. 

It is important to emphasize here that the (prompt) 
neutrino signal is expected to be associated in both time 
and direction with the electromagnetic emission of the 
source, and that this association can not be affected by 
the presence of inter-galactic magnetic fields, which may 
deflect the charged muons. The amplitude of a large 
scale inter-galactic magnetic field is limited by Faraday 
rotation measures to less than 1 uG [23|. Given the 
muon rest-frame lifetime, t^ = 2 x 10~^ s, the angle by 
which the muon may be deflected before decaying is only 
9 w tf,eB/m^,c ^ 10~^°{B/lnG), and the related time 
delay is Atg « 0^l^/c - l{B/lnG) (1^/10'^'^ cm) ms. On 
the other hand, neutrinos produced by the decay of sec- 
ondary muons (muons produced by high energy photons 



generated by inverse-Compton scattering of CMB pho- 
tons by e* pairs produced by the original high energy 
photon), may suffer strong delays since the mean free 
path of secondary electrons is large, ^ 10^^ cm, and they 
may therefore be significantly deflected by inter-Galactic 
fields. 

The comparison of arrival times of GRB-induced neu- 
trinos and GRB photons may allow one to test for 
Lorentz invariance violation (LIV) as expected due to 
quantum gravity effects [s], l2j]. Since the time delay 
is expected to be an increasing function of neutrino en- 
ergy, AtLiv « (rf/c)(e^/^£;pianck)", the IQi^ eV neutri- 
nos discussed in this paper would provide much more 
stringent tests than the pion decay ^^ 10^^ eV neutrinos 
expected to be produced within the GRB source. The 
accuracy with which the time delay could be measured 
is limited by the duration of the burst photons, which 
ranges from few ms's for short bursts to tens of seconds 
for long ones. A possible complication in this context is 
the possible emission of ^ 10^^ eV neutrinos during the 
GRB "afterglow" phase, which may take place on time 
scale of hours 25-|. For 10^^ eV neutrinos originating 
from z — \ [d ^ 10^^ cm), a time delay of hours would 
probe ^1 > 10* for n = 1 or ^2 ^ 10~^ for n = 2. It is 
important to note here that since the background level 
of neutrino telescopes is very low, especially at high neu- 
trino energy, the detection of a single EeV neutrino from 
the direction of a GRB on a time sale of even years af- 
ter the burst would imply an association of the neutrino 
with the burst and will therefore establish a time delay 
measurement. 
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